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Abstract 


A simple model to characterize the effects of heavy charged particles on 
metal-oxide semiconductor field-effect transistors (MOSFET) is presented. The 
model is applied to experimental results provided by the Air Force Weapons 
Laboratory, and an attempt made to simulate saturation phenomena observed in 
the threshold voltage change. The model assumes all trapped holes are within a 
few nanometers of the oxide-semiconductor interface, and takes into account the 
resultant counter electric field, and its effect on the yield fraction escaping recom¬ 
bination. An equation relating threshold voltage change as a function of dose is 
derived and used, but does not duplicate the saturation effect. This is because 
charge trapped near the interface reduces the internal oxide electric field very lit¬ 
tle. However, similarities in p-channel and n-channel MOSFET damage sensitivi¬ 
ties indicate considerable bulk charge trapping. This would modulate the internal 
electric field considerably. To determine whether this accounts for saturation will 
necessitate refinement of the model, taking hole transport and bulk trapping into 


account. 














MODELING THE EFFECTS OF HEAVY 
CHARGED PARTICLES ON METAL-OXIDE- 
SEMICONDUCTOR FIELD EFFECT TRANSISTORS 

I. Introduction and Background 

One of the most important devices used in very-large-scale integrated circuits 
is the metal-oxide-semiconductor field-effect transistor (MOSFET). Figure 1 
shows a schematic of an n-channel MOSFET. The device works as follows. 

When a positive voltage is applied to the gate, a strong electric field sets up 
through the gate oxide and penetrates into the p-type Si substrate. In the region 
just beneath the gate oxide, this field repels the majority carriers (which, for n- 
channel MOSFETs, are holes) and attracts minority carriers (electrons), thus 
forming an inversion layer. With a potential difference applied between source 
and drain, current flows through this inversion layer and the device is considered 
"on". The operation of a p-channel MOSFET is similar, only the p-type and n- 
type materials are reversed, and a negative voltage is applied to the gate contact. 

MOSFETs are very sensitive to the effects of ionizing radiation. The 
sequence of events associated with irradiation of an n-channel MOSFET is shown 
in Figure 2. Ionizing radiation, such as gamma-rays, electrons, and heavy charged 
particles, initially create many electron-hole pairs throughout the gate oxide. A 
certain fraction of these pairs are then annihilated through recombination. The 
positive applied voltage rapidly sweeps the remaining electrons from the oxide and 
into the metal gate contact. The holes transport more slowly toward the gate- 
semiconductor interface, where a certain fraction are caught in deep hole traps 
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Figure 1. Schematic of an N-channel MOSFET (Ref 1:22) 


and the rest disappear into the semiconductor. These trapped holes permanently 
alter the electric parameters of the device, primarily by lowering the threshold 
gate voltage required to create the inversion layer. In fact, if enough holes are 
trapped at the interface, an inversion layer may form without applying any vol¬ 
tage at all. For p-channel MOSFETs, the process is similar, except the holes 
migrate toward the metal-gate contact and are trapped at that interface. 

Recent experiments have been conducted by the Air Force Weapons Labora¬ 
tory to characterize the effects of ionizing radiation on MOSFETs. References 3 
and 4 describe this work in detail. Large numbers of non-radiation-hardened 
Intersil 3N161 p-channel and 3N171 n-channel discrete MOSFETs were irradiated 
with protons, electrons, and Co-60 gamma rays. After exposure to specified levels 
of radiation (under various bias conditions), the transistor gate threshold voltages 
were measured "in situ", and the changes in threshold voltage found as a function 
of radiation angle of incidence (the angle between the proton track and the electric 
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Figure 2. Sequence of Events Associated with Irradiation of MOSFET (Ref 
2:67) 


field), applied field, particle energy, and total ionizing dose. 

Of particular interest were the results of the proton irradiations. In this case, 
the MOSFETs were irradiated at room temperature with protons of energies from 
2-16 MeV. The angle between the proton track and the electric field was varied, 
and data collected for the cases of 0, -45, and 80 degrees. Figures 3 and 4 show 
the damage sensitivity (change in threshold voltage per dose) for the p-channel 
and n-channel MOSFETs, as a function of particle energy and angle of incidence. 
In general, the damage increases with both particle energy and angle of incidence. 
This is as expected, as will be explained in the theory section of this paper. 

What is not expected, however, is the anomaly where the damage caused by 


the 45 degree protons is greater than that of the 80 degree protons for particle 
energies less than 6 MeV. When first observed, the researchers thought that this 














TEST DEVICE 


Vf•: -S Volt* V4*i 0 of -10 Volt* 



PROTON ENERGY • M*V 


Figure 3. P-channel MOSFET Damage Sensitivity as a Function of Particle 
Energy and Angle of Incidence (Ref 4:1209) 

anomaly was due to the additional material (silicon dioxide passivation and 
aluminum metallization) the 80 degree protons had to pass through before reach¬ 
ing the gate oxide. If these additional material lengths were greater than the 
mean proton ranges, then the number of 80 degree protons reaching the gate oxide 
would be reduced. However, the actual distance traveled by the protons was 
found to be 6 to 7 times less than the mean ranges of 2-6 MeY protons. This 
implied that all of the protons reached and passed through the gate oxide, regard¬ 
less of the angle of incidence. 

In a further attempt to understand this anomaly, the data was replotted in a 
different format, shown in Figures 5, 6, and 7. Figure 5 shows the change in 
threshold voltage as a function of dose for the various angles of incidence. For 
this case, the gate bias potential was -5 Volts. Figures 6 and 7 show essentially 
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Figure 4. N-channel MOSFET Damage Sensitivity as a Function of Particle 
Energy and Angle of Incidence (Ref 4:1209) 


the same thing, though with different gate bias potentials. On the basis of these 
graphs, the researchers deduced that a saturation effect was occuring within the 
test sample, that had the same effect as reducing the internal electric field during 
irradiation. 

To explain this phenomena, the researchers proposed that the buildup of 
trapped holes at the interface resulted in a counter electric field. The yield frac¬ 
tion of holes escaping initial recombination increases with both the angle of 
incidence and the magnitude of the electric field. Initially, the higher yield frac¬ 
tion of the 80 degree protons results in a larger charge buildup at the interface. 
However, this larger buildup offsets the electric field within the oxide, decreasing 
the fraction of holes escaping recombination. Thus, the threshold voltage change 
due to 80 degree protons saturates before the 45 degree case. 
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Figure 5. Change in Threshold Voltage as a Function of Dose and Angle of Ir 
cidence for Applied Gate Voltage of -5 Volts (Ref 4:1212) 

The goal of this thesis project was to give this rough qualitative model a 
more precise quantitative formulation, and verify that the counterfield effect actu¬ 
ally explains the anomaly. The next section describes the stages shown in Figure 
2 in more detail. Assumptions used in the model are stated, and mathematical 
formula presented. The following section presents the major results and analysis, 
and is, in turn, followed by the conclusion. 
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Change in Threshold Voltage as a Function of Dose for 80 Degree 
Angle of Incidence and Gate Voltages of -5 and -10 Volts (Ref -4:1213) 
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II. Theory 


This section describes the assumptions and mathematical equations that 
enter into each stage of the model. Considered in turn are columnar phenomena 
including initial electron-hole formation and recombination, bulk phenomena 
including electron and hole transport, and interface phenomena including deep- 
hole trapping, and radiation-induced traps. Following this, an equation relating 
the change in threshold voltage to the total ionizing dose is developed. Both the 
dependence of the electric field on trapped charge, and the dependence of yield 
fraction escaping recombination on electric field are considered in this develop¬ 
ment. 


Columnar Phenomena 

As the protons pass through the oxide, they lose energy through ionization 
and excitation of electrons across the energy bandgap; electron-hole pairs are 
formed. The amount of energy needed to create an electron-hole pair in silicon 
dioxide has been estimated to be about 18 eV (Ref 5:1520). Protons, as well as 
other heavy charged particles, leave a very dense column of electron-hole pairs in 
their track. However, to determine the initial charge produced by an incremental 
dose, the electron-hole pairs are assumed to be evenly distributed throughout the 
bulk of the oxide. This incremental charge is given as 

AC=K g ADAlf v (1 

where 


AC = charge (Coulomb) 


















Kg = generation constant (Coulomb/cm 3 —rad) 

AD = dose (rad—Si0 2 ) 

A = area of oxide (cm 2 ) 

1 = length of oxide (cm) 
f y = fraction yield 

The generation constant, Kg, is found from 

_ (.01 J/kg—rad)(1.609 X 10 -19 C/hole)p 
^ (1.609 X 10' 19 )(1000 g/kg)W 

where 


p = density (g/cm 3 ) 

W = energy per electron—hole pair formed (eV/hole) 

For Si0 2 the density is 2.2 g/cm 3 and the energy required to form an electron hole 
pair is 18 eV. Filling these numbers into Eq (2) gives a generation constant of 
1.222 X 10~ 6 Coulomb/cm 3 —rad. 

Immediately after the electron-hole pairs are produced in the dense column, a 
certain fraction of them will recombine. The fraction that survive. f v . is a func¬ 
tion of the particle type and energy, the electric field within the oxide, and the 
incident angle between the particle track and electric field. To find the fraction 
yield, a numerical procedure developed by Oldham (described in references 6 and 
7) is used. A description of this procedure follows. 

Figure 8 shows a schematic of the track of a charged particle after passing 
through the silicon dioxide insulator. Looking down the direction of the track, 
the initial charge density configuration is assumed to be 


n 


± 



(3) 


where 
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Figure 8. Charged Particle Track Through the Oxide 
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n ± = electron—hole density (cm -3 ) 

N 0 = initial linear density (cm -1 ) 
b = Gaussian radius (cm) 
r = Distance center of the track (cm) 

The electrons and holes then move according to the equation (Ref 7:2695) 
dn ± dn ± 

-fa- = D ±^ n ± + “ Qn + n - 

where 

D ± = Diffusion coefficient (cm' s) 

= mobility (cm 2/ V—s) 

E x = Component of electric field in the x—direction (V/cm) 
a = Recombination coefficient (cm 3/ s) 

The first term on the right hand of Eq (-1) represents diffusion, the second term 
represents drift under the influence of the external field, and the third term 
represents recombination. 
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The actual program (found in the Appendix) sets up a grid in the x and in¬ 
direction with z in the direction of the proton track, and y perpendicular to the 
electric field for all proton angles of incidence. The initial electron-hole density is 
found from Eq (3) for each point in the grid. These densities are then placed into 
a finite difference form of the right hand side of Eq (4), and dn ± /dt found at each 
point on the grid. The new densities are then found from 

<9n. 

M^new) = n ±( r ^old) + ^~ At ( 5 ) 


where At is the increment in time. These new densities are inserted into the right 
hand side of Eq (4), and the process repeated until the electron and hole columns 
separate under the influence of the external field. Integrating n + (r,t) over the area 
of the grid gives a new linear hole density N(t). 

Figure 9 shows the electron-hole columns moving past each other under the 
influence of an electric field of arbitrary direction. For the case of a particle track 
perpendicular to the electric field, the fraction yield, f y , will be given simply by 
N(t)/N 0 , where N(t) is the linear density after the columns have separated. The 
analysis up to this point is only good for the region between the dashed lines in 
Figure 9, but it can be generalized to three dimensions. If it is assumed that the 
hole column to the right of the right-hand dashed line remains unchanged, then 
the fraction yield can be found from 


T 

/N(t)dt 


fy N T 


( 6 ) 


where T is the time required for the right hand of the electron column to cross the 
left hand of the hole column. 


The parameters /j ± , D ± , a, N 0 , and b in Eqs (3) and (4) are found in the fol¬ 
lowing way. The room temperature mobilities of silicon dioxide are found quite 


















Figure 9. Electron and Hole Columns Moving Under the Influence of Arbitrary’ 
Electric Field (Ref 7:2696) 

readily from the literature. The values used in this analysis were 
M + =2X 10 _s cm 2/ V—sec (Ref 8:2012) and - 20 cm 2/ Y-sec (Ref 9:1333). The 
diffusion coefficients, D ± , are then found from the mobilities by using the Einstein 
relation (Ref 7:2696) 



where k is the Boltzmann constant. T is the temperature, and e is the electron 
charge. The recombination coefficient, a, is found by using an expression derived 
by Langevin (Ref 7:2696) 

(A*+ + //_)e 

Q = - (8 

te 0 

where ( is the dielectric constant of silicon dioxide, and e 0 is the permittivity of 
free space. The initial linear density, N 0 , depends on particle energy and can be 
found from stopping power tables (Ref 10:310). The value of the Gaussian radius, 
b, was taken from Reference 7 where it was used as a variable meter to tit 


















this recombination model to experimental data. The value was found to be 3.5 
nanometers (Ref 7:2697). 

It should be pointed out that this model is only valid for situations where the 
electron-hole pair separation distance is short enough that the pairs lose their 
identity. In this case electrons readily recombine with holes from different pairs. 
This is generally true of the particle tracks left by heavy charged particles, 
because of their high linear energy transfer. For electrons and gamma-rays, the 
pair separation distance is large enough that recombination is pretty much res¬ 
tricted to electron-holes of the same pair. The model developed to handle these 
cases is called the geminate model (Ref 1:27). 

In general, the yield fraction f depends on three factors. First, it depends on 
the energy of the particle. This is because the initial linear density, N 0 in Eq (3). 
generally decreases with increasing particle energy, at least in the ranges of 
interest here. The lower the initial density, the less probability there is that 
recombination will occur, thus the yield fraction increases. Second, yield fraction 
also depends on the strength of the electric field. The stronger the field, the more 
quickly the electron and hole columns separate, decreasing the time available for 
recombination. Thus, yield fraction increases with increasing electric field. Third, 
the yield fraction depends on the angle of incidence. This can be understood by 
considering the two extremes. When the particle track is perpendicular to the 
field, the two columns are separated rapidly, which increases yield fraction. When 
the particle track is parallel to the field, the columns are pulled through each 
other, enhancing recombination, and reducing the yield fraction. Thus, yield frac¬ 
tion increases with increasing angle of incidence. 

This increase in yield fraction with particle energy, electric field and angle of 
incidence is the cause of the general characteristics of the curves seen in Figures 3 
through 7. In these figures, the change in threshold voltage also increases with 
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particle energy, electric field, and (except for particle energies below 6 MeV) angle 
of incidence. As will be seen below, the change in threshold voltage is directly 
proportional to the charge accumulated at the interface of the oxide and thus the 
yield fraction. 

Bulk Phenomena 

After escaping initial recombination the remaining electrons are rapidly swept 
from the oxide by the applied gate field. This process takes a few picoseconds to 
complete. Two assumptions are made about this stage. First, no electrons are 
trapped within the oxide. Second, having escaped columnar recombination, the 
electrons do not recombine further with holes previously trapped in the oxide. 

.After the electrons are swept from the device, the holes transport toward the 
interface much more gradually. .Assume all the holes move with the same velo¬ 
city, v = ^ + E. For an oxide of length 280 nanometers and an electric field of 
1 X 10 3 Y/cm at room temperature, it takes a hole 1.-40 X 10~ 5 seconds to cross 
the oxide. This time is short enough to justify ignoring any transient effects 
caused by transport. To simplify the analysis then, all the holes produced by an 
incremental dose, AD (that escape recombination), are assumed to instantaneously 
transport to the interface. Bulk trapping within the oxide is ignored. 

Interface Phenomena 

.As the holes pass through the interface region, a certain fraction, fj, are 
caught in deep hole traps within about 10 nanometers of the interface. The frac¬ 
tion trapped depends on the temperature and electric field, but is mostly a result 
of the processing techniques used during manufacture. The interface region has 











high local strain and a deficiency of oxygen atoms. This gives rise to a large 
number of Si-Si bonds instead of the usual Si-O-Si bonds found in silicon dioxide. 
When a hole passes one of these bonds, it can break the bond and recombine with 
one of the bonding electrons. This results in a Si atom retaining the remaining 
electron from the broken bond, and the positive charge residing with the other Si 
atom (Ref 1:34). 

Since the trapping fraction depends primarily on the manufacturing process, 
it can vary quite radically between different MOSFETs. Generally, f T can be as 
little as 1-2 percent for hard oxides, 10-20 percent for good quality commercial 
oxides, and as much as 50-70 percent in very soft commercial oxides (Ref 1:33). 

The spatial dependence of the charge trapped at the interface is also depen¬ 
dent on the relative hardness of the device. Hard oxides have distributions that 
fall off exponentially as one moves away from the interface into the oxide. On the 
other hand, the softer oxides have distributions that remain relatively constant 
(Ref 11:1207). 

These trapped holes can persist from a few hours to years. Annealing does 
take place, in which electrons from the Si semiconductor tunnel into the oxide and 
recombine with the trapped holes. However, since the time frame for annealing is 
considerably longer than the experiment, it is ignored in this model. 

This model assumes that a certain fraction, f T , of the holes escaping recombi¬ 
nation will be trapped in a constant charge distribution, p , extending a distance X 
from the interface. The incremental charge density due to an incremental dose is 
then 


v „ _ KgADIf y f x 


( 9 ) 


in the region near the interface and 0 etsewhere. X is presumably about 10 nm 


(Ref 1:33). The trapping fraction is uncertain, but since these devices are reported 
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as non-radiation hardened (Ref 3:4393), it is probably greater than 10 percent. 
This quantity is left undetermined in the analysis, but is assumed constant. 

Radiation-induced interface traps exist very close to the interface, and within 
the Si band gap. Their occupancy is determined by the position of the Fermi 
level at the interface, and can be negative (usually with n-channel devices), posi¬ 
tive (usually with p-channel devices), or neutral. In either case, the presence of 
these trapped charges results in shifts to the threshold voltage. At the present 
time, the precise mechanisms responsible for the formation of these states are not 
well understood, and a subject of considerable debate (Ref 1:37). These states are 
difficult to model. Since the deep hole traps are dominant (especially at the lower 
dose levels applicable here), interface traps are ignored in this analysis. 


Threshold Voltage vs Dose 

In this section, an analytical expression relating the change in threshold vol¬ 
tage to the total ionizing dose is found, so that the curves in Figures 5 through 7 
can be simulated. This expression will include the dependence yield fraction has 
on electric field and incident angle. 

The change in threshold voltage for arbitrary charge density is given by (Ref 
12:199) 

AV fh = J^^-dx (10) 

o f 

where 

1 = the length of the oxide (crn) 
p(x) = density of holes (cm -3 ) 
c = permittivity of Si0 2 


Assuming the charge density is constant between I — X and 1. and zero elsewhere. 

















the threshold voltage change becomes 


AV lb = ^MI (., 

or 

d AV th - ^ - m (12 

tn 2e 

where dAV th is the incremental change in threshold voltage resulting from an 
incremental change in charge density. Substituting the value for Ap given by Eq 
(9) into Eq (10) gives 

dAV lh - Kgfyffl (X - ^ 21) AD (13) 

In Eq (13), everything is assumed constant except the yield fraction. f v , which 
depends on the electric field, which in turn depends on the charge distribution. 

We next derive an expression relating the internal electric field to p and thus 
AV th . Figure 10 shows a schematic cross section of a MOSFET, with an areal 
density Q x < located at arbitrary x' within the oxide. This charge induces a charge 
density Qo at the gate, and Qj at the semiconductor. Ignoring any difference in 
work function between the metal gate and semiconductor, it can be shown that 

1 - \ f 

Qo = --—p—Q X ' (u: 

and 

Qi = —~pQx (1«V 

Using Gauss’ Law, the internal electric field due to charge density Q x . becomes 

E(x) = _ 'Sr^ Q » (16 ’ 


for 0 < x < x' and 


















Gate 


Qo u 


Oxide 


Qx 


Semiconductor 


Q, 


Figure 10. Schematic Cross Section of MOSFET With Arbitrary Areal Charge 
Density Q x . at x' 


E(x) - ~-Q x - 

for x' < x < 1. 


(17) 


Eqs (16) and (17) can be generalized to an arbitrary charge distribution p(x) 
as shown in Figure 11. Taking into account the field due to the external gate vol¬ 
tage and letting Q x = p(x' )dx' , the electric field at x becomes 


E (x) = KrtilW -n' + It 
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Expanding the second term on the right and rearranging gives 
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Figure 11. Schematic Cross Section of MOSFET with Arbitrary Charge Densi¬ 
ty 


The first term on the right hand side of Eq (19) is merely — AV th /l. Concern¬ 
ing ourselves with the field only at x < 1 — X, the second term becomes —pX/e. 
Finding p in terms of AV th from Eq (11) and inserting into this second term gives 


E(x) 





(20 


For simplicity, the yield fraction is assumed linear with electric field and can 
be given the form 

f v (E) = a*E + b, (21 

where a* and b^ will depend on the angle of incidence for a proton of given energy 
Substituting Eq (20) into Eq (21) will give f v as a function of change in threshold 
voltage. Placing f y (AV th ) into Eq (14) yields an equation of the form 
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III. Results and Analysis 


The previous section developed the model used to simulate the effect charged 
particle radiation has on MOSFET devices. This section describes the results of 
applying this model to the actual experiment described in the introduction and 
background section. 

Figure 12 shows the yield fraction as a function of particle energy for angles 
of incidence of 0, 45, and 80 degrees using the columnar recombination model. An 
electric field of 1.79 X 10 s V/cm was used. This field was found by dividing the 
gate voltage (5 Volts) by the thickness of the oxide (280 nanometers). The fraction 
yield increases with both particle energy and angle of incidence, and the charac¬ 
teristics of these curves are similar to those shown in Figures 3 and 4 except for 
the data anomaly. 

In finding a functional relationship between the change in threshold voltage 
and the dose (accounting for oxide field modulation due to hole trapping near at 
the interface) it was assumed that the fraction yield was a linear function of the 
electric field within the oxide. Figure 13 shows the fraction yield as a function of 
electric field for 4 MeV protons with angles of incidence of 45 and 80 degrees. 

This figure was found by using the columnar recombination model. As can be 
seen, the fraction yield increases with both electric field and angle of incidence, 
but the curves are not linear throughout the range of interest. To properly use 
Eq (23). a linear approximation was taken from 0.1 MY/cm to () 05 MY/cm (the 
general range of the oxide electric field during the irradiation). The coefficients in 
Eq (21) are found to be a, = 4.11 X 10~ 7 cm/V and b # = .0539 for an angle of 45 
degrees, and a$ = 4.88 X 10 7 cm/V and b$ = .0650 for an angle of 80 degrees, 
b'sing these values in Eq (23) with X = 10 nm, I = 280 nm. and V = 5 Volts, 
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Figure 12. Fraction Yield as a Function of Particle Energy and Angle of In¬ 
cidence. 

AV th was found as a function of dose for trapping fractions of .20. and .10. The 
results are shown in Figure 13. .As expected, the threshold voltage rises with dose, 
angle of incidence, and trapping fraction. However, the sought after saturation 
effect is missing. 

The reason the saturation effect is absent is straightforward. If all of the 
charge is assumed to accumulate within 10 nm of the interface, the charge giving 
rise to a threshold voltage change equivalent to the gate voltage will (from Eq 
(20)) only decrease the internal field to 0.982 of its original value. This is gen¬ 
erally true for charges at the interface—they modulate the internal electric field 
very little, certainly not enough to cause tlm duration effect seen in Figures '> 
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Figure 13. Fraction Yield as a Function of Electric Field and Angle of Incidence 

through 7. 

This raises the question whether all of the charge is actually trapped near the 
interface. There is evidence that it is not. Figures 3 and 4 show the damage sen¬ 
sitivity for both n-channel and p-channel MOSFETs. Under applied bias, holes 
created by radiation accumulate at the gate interface in p-channel MOSFETs. and 
at the semiconductor interface in n-channel MOSFETs. From Figures 3 and 4. it 
can be seen that the damage sensitivity is a little bit greater in the n-channel dev¬ 
ices. but not much. Our original assumption was that the charge was trapped 
within lOnm of the interface. Using Eq (10). a given charge, p, at the interface of 
the p-channel device results in a threshold voltage change of —50.p/e. If this same 
charge is distributed at the semiconductor interface of the n-channel device. Eq 
(10) gives a AV th = —2750p/f, which is 55 times greater. Of course this assumes 
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Figure 14. Threshold \'oltage Change As a Function of Dose. Angle of Incidence, 
and Trapping Fraction 


the trapping fractions are equal at both interfaces, which is unlikely. However, it 
seems improbable that the trapping fractions are that different. If a charge p is 
evenly distributed throughout the bulk of the oxide, then the change in threshold 
voltage is equivalent for both types of MOSFETs. 

Thus it appears that there is a significant amount of charge trapped in the 
bulk of the oxide. The hole transport process is actually very dispersive. Some 
holes move out of the oxide quickly, while others take a long time to do so. The 
assumptions made about hole transport trapping are generally true for thin oxides 
less than 100 nm and high electric fields, greater than 1 MV/cm. For the oxides 
considered here, the thickness was 280 nm and the field 1.79 X HF \ /cm. With 
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these parameters, it is likely that a significant proportion of charge is trapped in 
the bulk of the oxide (Ref 13:3940). 

The implication of large bulk trapping can be seen from Figure 15. This 
figure shows the electric field as a function of distance for a charge density distri¬ 
buted evenly throughout the oxide. The charge density is assumed large enough 
to give rise to a threshold voltage change equal to the applied gate voltage. From 
Eqs (11) and (19) the internal electric field is given as 

E(x)-2V.-L (24) 

8 l 2 

so E(0) = 0 and E(l) = 2V_/I. The distributed charge causes greater modulation 
throughout the bulk of the oxide, than that charge accumulated at the surface. 

To the left of the midpoint in Figure 15 the electric field falls significantly, 
causing a reduction in yield fraction. On the other hand, the field rises 
significantly to the right of the midpoint, which increases the yield fraction. 

From Figure 13, the yield fraction seems to slope off with decreasing electric field, 
so the decline in fraction yield in the left hand region may dominate the increase 
in the right hand side, resulting in an overall reduction in fraction yield. Whether 
this reduction is great enough to cause the saturation effect remains questionable. 

It could well be a contributor, but there are other things to consider. 

For instance, it was assumed earlier that the electrons, having escaped colum¬ 
nar recombination, were swept from the oxide without further recombination. 

For low hole densities and high electric fields this assumption is certainly valid. 
However, the cross section for electron-hole recombination is proportional to E - ’ 2 
(Ref 14:3203). The reduced field in Figure 15 could result in significant recombi¬ 
nation in the region to the left of midpoint. This could be another source of the 
saturation effect, especially if the field reaches zero within the oxide. Then the 
electrons would stop and recombine with a hole (Ref 16:1522-1523). 








Figure 15. Electric Field Within the Oxide from an Evenly Distributed Charge 
Density 

A third assumption made in this mode! was that the fraction of holes 
trapped at the interface was constant. However, it is known that the hole trap¬ 
ping cross section is dependent on the electric field and varies with E -1 * (Ref 
15:1192). The increasing field at the interface, as seen in Figure 15 may result in 
a decline in trapping fraction, that could also be a contributor to the saturation 
effect. 

These three field modulated mechanisms—yield fraction, electron recombina¬ 
tion, and hole trapping may all contribute to the saturation effects seen in Figures 
5-7. To verify this will require a more refined model that not only takes these 
mechanisms into account, but also accounts for the wide lispersion inherent in the 
hole transport. This will be required to characterize the hole buildup in the bulk 
of the oxide. 


27 

















fV. Conclusion 


In conclusion, this simple model has shown that the electric field reduction 
due to holes trapped at the interface cannot cause the saturation effect. Charges 
trapped close to the interface cannot reduce the internal electric field to any great 
degree--certainly not to the extent required to effect the yield fraction of holes 
escaping columnar recombination. However, similar p-channel and n-channel dam¬ 
age sensitivities indicate that considerable charge is trapped within the bulk of the 
oxide. Charge distributed within the oxide bulk modulates the electric field con¬ 
siderably. This modulated electric field could contribute to the saturation effect in 
three ways— reduction of the yield fraction throughout the oxide, increase in elec¬ 
tron recombination, reduction of hole trapping at the interface. 

To verify this, the present model would have to be refined in a number of 
ways. First, the dispersive nature of hole transport would have to be considered 
so that reasonable estimates of the trapped hole distribution (both at the interface 
and within the bulk of the oxide) can be made. Second, Eqs (10) and (19) would 
have to be modified to numerically find the threshold voltage and electric field for 
an arbitrary charge density. Third, the nonlinearity between yield fraction and 
electric field would have to be accounted for, so that overall yield fraction due to 
a widely varying internal electric field may be found. Fourth, the electron tran¬ 
sport and electric field dependent recombination cross section would have to be 
introduced. Finally, the hole capture cross section dependence on electric field 
would have to be taken into account. Introducing these refinements into the 
model should allow determination of the primary cause of the saturation effect 
and explain the data anomaly. 
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Appendix: Columnar Recombination Computer Code 


This appendix contains the code used to find the yield fraction of electrons 
and holes escaping inita.il recombination. The following algorithm is used: 


Step 1. 


Step 2. 


Step 3. 


Step -4. 


Step 5. 
Step 6. 


Step 7. 
Step 8. 
Step 9. 
Step 10. 
Step 11. 

Step 12. 
Step 13. 


Input the initial linear density, the mobilities, the diffusion coefficients, 
the recombination coefficient, the electric field, the angle of inci 
incidence, the column radius, and the length of the oxide (lines 22-23). 

Set up the x and y grid. In this case the maximum values in the x and 
y directions and the number of grid spaces are read in. The program is 
set up to take advantage or the symmetry in the y direction (lines 30- 
45 ). 


Input and initialize time related elements. These include the time incre¬ 
ment and the ratio of time increments between each step. This ratio is 
inserted to help shorten the run times. The time variable is initialized 
and the maximum time calculated (lines 51-53). 

Set up boundary conditions on the grid. In this routine the boundaries 
are assumed large enough that the electrons and holes separate well 
before the columns reach the boundary time so the electron and hole 
densities at xrnax, xmin. and ymax are all set to zero (lines 58-73). 

Determine the linear density of particle at time equal to zero (lines 108- 
124). 


Determine the yield fraction for each time step (lines 132-216). 

Determine new time and time ratio (lines 133-134). 

Set new boundary condition along ymin (lines 139-1 12) 

Determine the new density at each point on grid (lines 144-183). 

Find the diffusion term in Eq (2) (lines 151-160). 

Find the drift term in Eq (2) (lines 164-170). 

Find the recombination term in Eq (2) (line 17 1). 

find the time rate of change of the holes electrons (lines 179- 
1 s(l). 


Step 14. 


find the new charge density at the new time (lines 188-196). 





















Step 15. Find the number of charge carriers per unit length and thus the 

new yield (lines 200-212). 

Step 16. Print out the time,electron and hole density, and the yield (line 

213). 

Appendix of variables. 

alpha - recombination coefficient in cm3/seo 
b - Gaussian radius in nanometers 
bcm - Gaussian radius in centimeters 

check - variable used to determine whether the point is outside the column or not 
delx - the length of each element in the grid in the x- direction in nanometers 
dely - the length of each element in the grid in the y- direction in nanometers 
delt - the time increment in seconds 

derele(i.j) - time rate of change of electron charge density in /cm3*s 
derhol(i.j) - time rate of change of hole charge density in /cm3*s 
diffn - diffusion term in the electron density rate of change equation 
diffp - diffusion term in the hole density rate of change equation 
dndx - first derivative with respect to x of electron density 
dpdx - first derivative with respect to x of hole den density 
dplus - diffusion coefficient of holes in cm2/s 

driftn - drift term in the electron density rate of change equation 
driftp - drift term in the hole density rate of change equation 
dneg - the diffusion coefficient of the electrons in cm2/s 
d2ndx2 - second derivative with respect to x of the electron density 

d2ndy2 - second derivative with respect to y of the electron density 

d2pdx2 - second derivative with respect to x of the electron density 

d2pdy2 - second derivative with respect to y of the hole density 














e - the electric fiefield applied across the oxide in Volts/cm 
eleden(i,j) - the electron charge density in charges /cm3 

eledenx - the density of electrons along the x axis 

eledeny - the density of electrons along the v axis 

elenum - the number of electrons along the y axis in charge/cm 

elenuml - the number of electrons above the y axis in charge/cm 
elenumo - the total number of electrons equal to elenum+elenum 1 *2 
hx - thethe length of each element in cm in the x direction 
hy - the length of each element in cm in the y direction 
holden(i,j) - the hole charge density in charges /cm3 
holdenx - the hole density along the x axis in charges/cm3 
holdeny - holthe hole density along the y axis in charges/cm3 
holnum - the number of holes along the y axis in charges/cm 
holnuml - number of holes above the y axis in charges/cm 
holnumo - total number of holes equal to holnum+holnum 1 *2 
length - length of oxide in cm 

m - number of elements between upper and lower boundaries in the grid in x 
direction 

midx - midpoint of the grid in the x direction 
midy - midpoint of the grid in the y direction 
muneg - mobility of electrons in cm2/V*sec 
muplus - mobility of holes in cm2/V*sec 

n - number of elements between the upper and lower boundaries of the grid in th 
y direction 

no - initial line density in charges/cm 

recoin - recombination term in the density rate of change equation 

thetad - angle of incidence between the particle track and the electric field in 
degrees 

thetar - angle of incidence in radians 



















time - time for current step 

tmax - maximum time for program to run 

tration - ratio of delt from one iteration to that of the previous one 

x(i) - values of x for each grid point i in nano- meters 

xmax - maximum value of x in nanometers 

xmin - minimum value of x in nanometers 

y(i) - values of y for each grid point i in nano- meters 

yield - fraction yield for columnar recombination 

vmax - maximum value of y in nanometers 

ymin - minimum value of y in nanometers 










0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 


c This code is designed to find the concentration of holes and elec- 

c trons(in charges/cm**3) and the yield fraction of holes escaping 

c recombination aue to interaction of charged particles (alphas and 
c protons) with silicon dioxide. 

implicit none 

real pi,no,muneg,muplus,dneg,dplus,alpha,e,thetad,thetar 
real b,bcm,xmax,xmin,ymax,ymin,delx,dely.x(201),y(102) 
real delt,eleden(201,102),holden(201,102),check.yield 
real derele(201,102),derhol(201,102),hx,hy,tmax.time.length 
real elenum.holnum 

real elenumo,holnumo,eledenx(201,100).holdenx(201.100) 

real eledeny(l01,100),holdeny(101.100),elenum 1,hoinum 1 ,d2ndx2 

real d2ndy2,diffn,dndx,driftn.recom.tratio 

real d2pdx2.d2pdy2,diffp,driftp.dpdx 

integer n,m,midx,midy,i,j 

open( unit= 1,file=' rec.dat'.st at us='old'.READONLY) 


c This section initializes a number of the constants. 

read! l,*)no,muneg,muplus,dneg, dpi us, alpha 
read(l,* je.thetad.b, length 
pi=3.1416 

thetar=thetad*pi/l80 


c This section sets up the x and y grid. 

read(l,*)xmax,xmin,ymax,ymin,m.n 

bcm=b*lE-07 

midx=m/2+l 

midv=2 

delx=(xmax-xmin)/ m 
dely=(ymax-ymin)/n 
hx=delx* IE-07 
hy=dely*lE-07 
x(l )=xmin 
y(l)=ymin 
do 10 i=2.m+l 
x(i)=x(i-l)+delx 
10 continue 

do 20 j=2,n+l 
y(j)=y(j-l)+dely 
20 continue 


c This section inputs the time related elements such as the time 
c increment, and the maximum time considered. 

read(l,*)delt,tratio 

tmax=(length/cos(thetar))/(muneg*cos(thetar)*e) 

time=0.0 
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c This section inputs the boundary conditions for the charge density 
c at xmin and xmax. 

do 30 j=2,n+l 
eledenfl,j)=0.0 
eleden(n+l,j)=0.0 
holden(l,j)=0.0 
holden(n+1 ,j )=0.0 
30 continue 


c This section inputs the boundary conditions for the charge density 
c at ymax. 

do 40 i=2,m 

eleden(i,n+l)=0.0 
holden(l,n+l)=0.0 
40 continue 


c This section initializes the charge densities at time equal zero. 

do 50 i=2,m 
do 60 j=2,n 

check=-(x(i)**2+y(j)**2)/b**2 
if(check.lt.-100.00)tnen 
eleden(i,j)=0.0 
holden(i ,j)=0.0 
else 

eleden(i,j)=(no/(pi*bem**2))*exp(-(x(i)**2+y(j)**2)/b**2) 

holden(i,j)=eleden(i,j) 

endif 

60 continue 
50 continue 


c This section finds the charge density along the x and y axis at 
c time equal to zero. 

do 70 i = 1, m+1 

eledenx(i.l)=eleden(i,midy) 

holdenx(i,l)=ho!den(i,midy) 

70 continue 

do 80 j = 1 ,n + l 

e!edeny(j,l)=eleden(midx.j) 

holdeny(j,l)=holden(midx.j) 

80 continue 


c This section determines the total number of charge carriers(per 
c unit length) at time equal to zero. 

elenum=0.0 
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holnum=0.0 
elenuml=0.0 
holnuml=0.0 
do 90 i=2,m 
do 100 j=3,n 

elenuml=elenuml+eleden(i,j)*(hx*hy) 
holnuml = holnumH-holden(i,j)*(hx*hy) 
100 continue 
90 continue 

do 110 i=2,m 

elenum=elenum+eleden(i,2)*hx*hy 
holnum=holnum+holden(i,2)*hx*hy 
110 continue 

elenumo=elenum-l-elenuml*2 
holnumo=holnum+holnuml*2 
yield—holnumo* del t 
print*,time,elenumo,holnumo,yield 
print*,'max time=',tmax 


c This section determines the time derivative of the charge density 
c the new charge density and the total charge as a function of time. 

120 if(time.gt.tmax)go to 130 
time=time+delt 
delt=delt*tratio 
elenum=0.0 
holnum=0.0 
eienuml=0.0 
holnuml— 0.0 
do 140 i = l,m + l 

eleden(i,l)=eleden(i,3) 

holden(i,l)=eleden(i,3) 

140 continue 

do 150 i=2,m 
do 160 j=2,n 


c This section finds the diffusion term for the electrons from finite 
c difference form. 

d2ndx2=(eleden(i+l,j )-2*eleden(i,j)+eleden(i-1 ,j ))/hx* *2 
d2ndy2=(eleden(i,j + I )-2*eleden(i,j )+eleden(i.j-1 ))/hy**2 
diffn=dneg*(d2ndx2+d2ndy2) 

c This section finds the diffusion term for the holes from finite 
c difference form. 

d2pdx2=(holden(i + l ,j l-2*holden(i,jf+holden(i-1 ,j))/hx** 
d2pdy2=( holden(i.j +1 }-2* holdenf i.j )-*-h<'lden( i.j-1 ))/hv* * 
diffp = dplus*(d2pdx2 + d2pd}2) 

r This section finds the drift term for the elect i 
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ci ci 















dndx=(eleden(i + l,j)-eledenfi-l,j))/(hx*2) 
driftn=muneg*e*sin(thetar)*dndx 

This section finds the drift term for the holes. 

dpdx=(holden(i + l,j)-hoIden(i-l,j))/(hx*2) 
driftp = muplus*e*sin(thetar)*dpdx 

This section finds the recombination term for both holes and electrons 

recom=-alpha*eleden(i,j)*holden(i,j) 

This section then finds the time rate of change of the electron and 
hole density for each x and y. 

derele(i,j)=diffn+driftn+recom 

derhol(i,j)=diffp+driftp+recom 


150 continue 

c This section finds the new charge density at each point for 
c new time. 

do 170 i=2,m 
do 180 j=2.n 

eleden(i,j)=eleden(i.j)-Hlerele(i,j)*de)t 

holden(i.j)=holden(i.j)+derhol(i,j)*delt 

if(eleden(i,j).lt.0.0)then 

eleden(i,j)=0.0 

endif 

180 continue 

170 continue 

c This section finds the number of charge carriers(per unit length). 

do 190 i=2,m 
do 200 j=3.n 

elenum 1 =elenum 1 +eleden(i.j)*(hx*hy) 
hoi num 1 = hoi n urn 1 +holden(i ,j )* (hx * hy) 

200 continue 

190 continue 

do 210 i=2.m 

elenum=elenum+eleden(i,2)*hx*hy 
holnum = holnum + holden(i.2)*hx*hy 
210 continue 

elenum=elenuml*2+elenum 
holnum = holnum.l *2+holnum 
yield=vield+holnum*delt 
print*, time,elenum. hoi num,yield 

go to 120 
130 continue 
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19. Abstract 

A simple model to characterize the effects of heavy charged particles 
on metal-oxide semiconductor field-effect transistors (MOSFET) is presented. 
The model is applied to experimental results provided by the Air Force Wea¬ 
pons Laboratory, and an attempt made to simulate saturation phenomena ob¬ 
served in the threshold voltage change. The model assumes all trapped holes 
are within a few nanometers of the oxide-semiconductor interface, and takes 
into account the resultant counter electric field, and its effect on the frac¬ 
tion escaping recombination. An equation relating threshold voltage change as 
a function of dose is derived and used, but does not duplicate the saturation 
effect. This is because the charge trapped near the interface reduces the in¬ 
ternal oxide electric field very little. However, similarities in p-channel 
and n-channel MOSFET damage sensitivities indicate considerable bulk charge 
trapping. This would modulate the internal electric field considerably. To 
determine whether this accounts for saturation will necessitate refinement 
of the model, taking hole transport and bulk trapping into account. 
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